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Abstract: Using the example of selected decays driven by R-parity- violating super- 
symmetric operators, we demonstrate how strong final-state interactions can be controlled 
quantitatively with high precision, thus allowing for a more accurate extraction of effec- 
tive parameters from data. In our examples we focus on the lepton-flavor-violating decays 
r — > fiir + ir~ . In R-parity violation these can arise due to the product of two couplings. We 
find bounds that are an order of magnitude stronger than previous ones. 
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1 Introduction 

If weak-scale supersymmetry [1] is the solution to the hierarchy problem there must be new 
supersymmetric particles with masses below 0(5 TeV), which are accessible to the LHC. 
There are basically two forms of the minimal supersymmetric Standard Model, with light 
supersymmetric fields, distinguished by their superpotential. The most widely studied case 
is R-parity conservation, where the symmetries of the supersymmetric Standard Model are 
extended to include the discrete multiplicative symmetry R-parity. This renders the proton 
stable in the theory and the resulting renormalizable superpotential is given by 



WMSSM = tab 



(h e )ijLiHjEj + (h d )ijQiH^Dj + (h u )ijQ1H u Uj + nH%H b 



(1.1) 



Here L, Q,Ha, H u are the lepton, quark, and Higgs SU(2)x doublet left chiral superfields, 
respectively. E, D, U are the corresponding SU(2)x singlet lepton and quark left chiral 
superfields. i,j,k £ {1,2,3} are generation indices and a, b E {1,2} are SU(2)l gauge 
indices. The hf= e du are dimensionless Yukawa couplings and [i is the Higgs mass mixing 
term. Note that if non-renormalizable terms are allowed then for example the R-parity 
conserving superpotential term QQQL leads to a dimension-five proton decay operator [1]. 
This is thus only suppressed by one power of the large mass scale. In this case proton 
hexality [2, 3] is the appropriate symmetry. It leads to the same low-energy superpotential 
given in Eq. (1.1), but prohibits all dimension-five proton decay operators. 

Alternatively, the proton is also stable for the discrete Z3 symmetry baryon triality [4]. 
The resulting renormalizable superpotential is [5] 



Wb 3 = Wmssm + e, 



ab 



\jkLiL b jE k + \' ijk LfQ b j D k + KiLlH b a 



X2) 
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Here the A^, A- jfc are dimensionless couplings and the have dimension mass. Both 
baryon triality and proton hexality are discrete gauge anomaly- free in the sense of Refs. [2, 
4, 6]. At any given energy scale the Ki can be rotated to zero, by a transformation in 
(Li, Hd) space [7, 8]. Since we perform our computations at a fixed low-energy scale, we 
shall focus on the case Ki = in the following. 

An important feature of the additional terms in Eq. (1.2) is that they violate lepton 
number and flavor. Correspondingly there is a large set of bounds on the couplings Xijk 
and A^- fc . This was first considered in Ref. [7], including also possible contributions to 
rare meson decays. A more systematic approach was taken in Ref. [9], considering bounds 
on all couplings individually. Since then many bounds have been set from the decays of 
mesons, including also bounds on products of operators [10-18]. In Ref. [9], for example, 
the authors considered the R-parity-violating contributions to 

t^w^v < L3 > 

1 (7T — > \XV) 

An operator \' 2lk L2Q\D] i contributes only to r(7r — > fiv), thus modifying R w . In computing 
the exact contribution to r(-7r — > fj,u), one uses the definition 

(0\u-fP L d\ir-(q)) = ~f n q» (1.4) 

of the pion decay constant f n for the V — A current, and the current-algebra approximation 

(0|uP L d|7rfa)> = V2i M * f n (1.5) 
m u + m d 

for the chiral (pseudo) scalar coupling. Here, M n denotes the charged pion mass, m u d are 
the first-generation quark masses, and we make use of the left- and right-handed projection 
operators Pl/r = (1 T 7s)/2- Using the experimental value including the error of R n [19], 
which agrees with the Standard Model, results in a bound on A^. This assumes that 
LiQ\Dk is the sole new operator contributing. 

More recently Herrero and collaborators have studied r decays in various R-parity 
conserving supersymmetric models [20, 21], for example the decays 

r ± -»• ^ ± + / /vr + ^~, it it /K + K~, K°i? /^, v' / P°, • (1-6) 

The authors go beyond the simple current algebra approximation and employ chiral per- 
turbation theory and resonance chiral theory [22]. They thus dramatically improve the 
precision of the computation and therefore also the resulting bounds on new physics. The 
purpose of this paper is to refine these techniques further — in particular, we include both 
the scalar form factors and the vector form factor, model-independently. In addition, we 
will also discuss the case of R-parity violation. Specifically we shall focus on the decay 

-> /^tt+tT (1.7) 

to clarify our method. In terms of R-parity-violating operators, this decay receives contri- 
butions via the parton level processes 

— > pL UiUf- and — > [i^djd^ . (1.8) 
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Figure 1. Fcynman diagrams inducing transitions of the type t~ — > fi~uu (a) or t~ — > fi~dd 
(b)-(d) by exchange of supersymmetric particles. The arrows denote the direction of the flow of 
the left-handed fields. The particle labels denote the incoming (t — ) or outgoing fields (all others). 
Thus in (b), e.g., we have an outgoing (SU(2) singlet) dj and (SU(2) singlet) di . 

Combining the operators e a bXijkL'}L b jEf t and eabX^^L^Q^Dk from the superpotential in 
Eq. (1.2) and integrating out the heavy intermediate scalar fermion, we obtain several 
independent contributions to the decay r~ — > pT + 2 quarks (cf. Appendix A): 



(a) Combining e a bX' 3ij L^Q^Dj and [e a bX' 2k .L^Q^D^ corresponds to r 



2fe r 

1 A' A'* 

= - 3lJ 2 2 3 (u k 'y a P L Ui)(fi>y a PLT) . 

2 m% 



\x UiUk with 
(1.9) 



(b) Combining e a bX' 3i jL'i i Q b i Dj and [e^A^Lf Q\Dg\ ^ corresponds to r 



eflf 



1 Kij^2U 

2 ml 



(d j j a P R d e )(ii>y a P L T) . 



(c) Combining e ah \yaL\L\Ei and \e ah \\- k L a i Q h -D } )\ i corresponds to r" 



VXmKjk {dtPRd^PLT) . 



2 mj 



► fi djde with 
(1.10) 

jj,~djdk with 
(1.11) 



(d) Combining 6^X2^2^^ and [e a6 A- ^L^Q^D^ corresponds to t ->• fi djdk with 



J?cs = \^^(d k P L d 3 )(jiP R T) . (1.12) 
2 m~ 

Diagrammatic representations of contributions (a)-(d) are shown in Fig. 1. 
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2 Application to the vector current 



The invariant mass distribution of the width for the decay r — > fj>ir + ir is given by 

dF 1 1 



where s is the invariant mass squared of the pion pair, |p*+| (^*+) is the momentum 
(angle) of the ir + in the rest frame of the pion pair, and |p^| and fi^ are to be given in the 
r rest frame. 

The essential observation is that the matrix element A4 factorizes, since the primary 
transition is short-ranged (the range of interaction is set by the inverse mass of the ex- 
changed super symmetric particles), while the final-state interaction is long-ranged. There- 
fore, we may write 

M = (7r + 7r-|dM|0>M { r Q) . (2.2) 

The reduced matrix elements are to be calculated in the underlying, fundamental 

theory, while the hadronic matrix elements (7r + 7T can be deduced either from 

data or determined with theoretical input. 

To begin with, let us assume that only the vector current contributes here — the gen- 
eralization to the scalar current is straightforward and will be presented below. As long 
as isospin is assumed to be conserved, two pions with vector quantum numbers (i.e., in 
a P-wave) only couple to the isovector component of the current. Therefore we need to 
consider only a single form factor and the hadronic matrix element is given by the pion 
vector form factor, Fy(s), defined via 

<7r+(p 7r+ )7r-(p 7r -)||(u 7 ^- drrd)\Q) = F v {a){p^ - P7T -) a . (2.3) 

Fy(s) is very well known both from direct measurements of e + e~ — > 7r + 7r~ [23-28] and, 
via an isospin rotation, of r~ — > ir~ir°v T [29], as well as theoretical studies [30-34]. It 
collects all non-perturbative irir interactions and is universal in the elastic region, which to 
excellent approximation comprises the energy range s < lGeV 2 . 1 

We briefly illustrate how to describe Fy(s) theoretically, based solely on the funda- 
mental principles of analyticity and unitarity. Figure 2 gives a graphical illustration of the 
discontinuity of the form factor, regarded as an analytic function of s in the complex plane 
cut along (parts of) the positive real axis, in the elastic regime, i.e., considering two-pion 
intermediate states only. It is given by 

disc F v (s) = F v (s + ie) - F v (s- ie) = 2i Im F v (s) 

= 2i F v (s) a{s) t\ (s)* 6(s- 4M 2 ) = 2i F v (s) sin 5{ (s)e~ i5 ' {s) 6{s - 4M%) , 

(2.4) 



1 In the following section, we shall emphasize the importance of the large inelastic coupling of the pion- 
pion isospin I — S'-wave to KK intermediate states in the region of the /o(980) resonance. We wish to 
point out here that the coupling to kaons has, in contrast, an entirely negligible effect on the pion-pion 
P-wave. In this case, the inelasticity is dominated by 4-7T intermediate states, often thought to be effectively 
clustered as ttlu (compare Ref. [35]), and only rises very slowly roughly above the ttlu threshold. 
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disc 



'X/X/X/^^ 



'X/X/X/^^ 



Figure 2. Graphical representation of the discontinuity relation for the pion form factor. The 
dashed lines arc pions, while the wavy line stands for the current coupling to these. The black disc 
denotes the form factor, while the gray disc denotes the pion-pion scattering T-matrix (which, by 
angular momentum conservation, is automatically projected onto the appropriate partial wave). 



where a(s) = \J\ — 4M% / s is proportional to the two-particle phase space, and tj(s) refers 
to the pion-pion partial-wave amplitude of isospin / and angular momentum £, obtained 
from expanding the corresponding T-matrix in Legendre polynomials. In the final step, 
we have rewritten the P-wave amplitude in terms of the phase shift <5{(s) in the canonical 
manner. One immediately deduces Watson's final-state theorem [36]: reality of Imi*V(s) 
implies that the phase of Fy(s) coincides with <5j(s). The solution to Eq. (2.4) is given by 

F v (s) = P v (s)n\(s) , n\(s) = exp\- [°° da' ^ ) , (2.5) 

where ^i(s) is the Omnes function [37] and Pv(s) a polynomial. The pion-pion phase 
shifts are known to excellent precision (up to at least y/s ~ 1.1 GeV) from analyses of the 
highly constrained system of dispersion relations known as Roy equations [38-41]. 

Perturbative QCD suggests Fy(s) should fall off like 1/s for large values of s (up to 
logarithmic corrections) [42], which is also the behavior of the Omnes function if the phase 
shift approaches it asymptotically, as phenomenology indeed suggests. Hence, -Py(s) is 
required to be a constant. Gauge invariance finally requires the normalization to be fixed 
to Fy(s = 0) = 1, therefore Pv{s) = 1. The representation Eq. (2.5) can be improved 
by taking inelastic effects (beyond two-pion intermediate states) into account. We here 
use a parametrization presented in Ref. [43], which describes the available high-precision 
data [23-29] perfectly. 

The relevant reduced matrix elements needed to complete Eq. (2.2) — read off from the 
effective Lagrangians given in Eqs. (1.9)-(1.12) — can be subsumed in the expression 

M Va = X V Mp»haPLT(p T )} , (2.6) 

where the effective coupling Ay is given by 

A 1/_^^_A^\ 
4 V m ~ mi, I 

This in principle comprises six different contributions = 1, 2, 3) which can enhance each 
other or lead to cancellations, depending on the phases of the R-parity-violating couplings. 
In the Standard Model there is a strong hierarchy among the Yukawa couplings. For 
example the top quark Yukawa coupling is almost a factor of forty larger than the bottom 
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[GeV] yfi [GeV] yfs [GeV] 

Figure 3. The predicted signals individually for the currents (a) qj a q, (b) (uu + dd)/2, and (c) ss. 
In all cases the effective coupling constant is set to 1 GeV -2 . For the uncertainty bands of the scalar 
form factors, see the discussion in the main text. 



quark Yukawa coupling. Since no supersymmetry with or without R-parity has yet been 
found, we shall for simplicity consider one product of operators at a time in Eq. (2.7). We 
thus employ the assumption that there is a hierarchy in the unknown R-parity-violating 
couplings [9, 44]. 2 

Using Eqs. (2.3) and (2.6) in the definition of the matrix element we find for the 
spin-averaged squared matrix element 



spins 



A V 



AM 2 ) 



2 i 2 

m T + m fl 



s + 



\M>)\' 



(2.8) 



where we included a prefactor of 1/2 from averaging the incoming r polarizations. 
X(x,y,z) = x 2 + y 2 + z 2 — 2{xy + xz + yz) is the usual Kallen function, and z = cos 8 
denotes the angle of the pions relative to the leptons in the 7r7r rest frame. Evaluating the 
angular integral and collecting all kinematic prefactors, we arrive at 



\2 



dTy _ 
d\fs 2567r 3 m^ 



(s-4M 2 f 2 X^ 2 (m 2 T ,ml,s) 



' A(m 2 ,m 2 ,s) 
"V+m^-sH 



\F v (s)\\ (2. 



Since all quantities but Ay are known in Eq. (2.9), it provides the kind of expression we 
are looking for. In particular, the norm of the hadronic current is fixed unambiguously in 
contrast to previous studies, where quark-model wave functions were employed to fix the 
normalization. The resulting y/s distribution is depicted in panel (a) of Fig. 3, where we 
have used Xy = 1 GeV -2 , which is easily rescaled. 

The signal in the vector channel can be represented reasonably well by a Breit-Wigner 
distribution with an energy-dependent width as provided by various parametrizations. Here 
the advantage of our approach lies mainly in its fixed normalization — we briefly compare to 



2 We keep in mind that Nature does not necessarily obey such analogies. For example the PMNS neutrino 
mixing angles are large compared to the CKM quark mixing angles. 
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the narrow-resonance approximation (where one studies the decay r — > // p (770), assuming 
a stable p ) in Appendix B. An additional advantage is that not only the modulus of the 
form factor, but also its phase is fixed unambiguously, such that the interference of various 
currents can be analyzed as well. 



3 Application to the scalar currents 

Contrary to the vector currents, scalar currents are typically not well represented by Breit- 
Wigner functions (for a detailed discussion see e.g. Ref. [45]). In this case, isospin symmetry 
requires two pions to couple exclusively to isoscalar scalar currents. However, there are 
two such form factors, originating from non-strange and strange isoscalar scalar sources. 
They are called and r* , respectively, and will contribute simultaneously. The use of 
scalar form factors (as opposed to Breit-Wigner parametrizations) is unavoidable if one 
wants to determine the underlying couplings in a controlled way. 
The expressions analogous to Eq. (2.3) now read 

{tv+tt- \qq\0) = B^Tlis) , (3.1) 

where the quark flavors may be either qq = {uu + dd)/2 for the light quarks, with the 
superscript q = n denoting the corresponding scalar form factor, or qq = ss for strange 
quarks (with superscript q = s). Furthermore, B n = M%/(m u + md), B s = {2M^ — 
M%)/(2m s ). With this convention, the form factors, T^(s), are invariant under the QCD 
renormalization group. For the numerical evaluation, we will use the values obtained from 
averaging lattice computations with Nf = 2 + 1 dynamical flavors [46], {m u + rrid)/2 = 
(3.43 ±0.11) MeV, m s = (94 ± 3) MeV, to be understood in the MS scheme at the running 
QCD scale /i = 2 GeV. In addition, in analogy to Eq. (2.6) we may now write 

MY = \l [fl(p,)P R T(p T )] , (3.2) 

where, again, the superscript q denotes the quark flavor fed by the corresponding operator. 
Comparison to Eqs. (1-11) and (1.12) yields 

\n _ (^3i2 + A 2 j 3 )A-i 1 _ (A 3 i 2 + A 2 j 3 )A-22 /o o\ 

Xs= 4m| ' Ks= 4m| " (3 ' 3) 

With these expressions we find 

^ = ^4^3 ( s ~ m D 1/2x1/2 K, ml, s) [(m T + m,f - s] \B n X n s T n n (s) + B s X s s K(s) | 2 . 

(3.4) 

Experimentally, the scalar form factors are not accessible as directly and unambigu- 
ously as the vector form factor. However, they can be reconstructed from dispersion theory, 
similar to what we discussed as the Omnes representation of the vector form factor in the 
previous section. The main difference is that the elastic approximation breaks down much 
earlier in the pion-pion 5-wave (of isospin 1 = 0) due to the strong inelastic coupling of 
two S*-wave pions to KK in the region of the /o(980), i.e., beginning immediately at the 
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KK threshold. In order to describe the scalar form factors including energies around the 
mass of the /o(980), it is therefore mandatory to solve a two-channel Muskhelishvili-Omnes 
problem [37, 47]. The discontinuity equation (2.4) has to be generalized to two coupled 
channels for the scalar form factors of pion and kaon. The two-channel T-matrix corre- 
spondingly can be parametrized in terms of three input functions, the pion-pion 5-wave 
phase shift 5q(s) known from Roy equation solutions [39-41], as well as modulus and phase 
of the inelastic reaction tttt — > KK [48-50]. We again assume the fall-off of all scalar form 
factors oc 1/s, and the scattering phases involved to approach integer multiples of n in 
the appropriate way. The solution to the coupled-channel discontinuity equation cannot 
be written down analytically in a similarly compact form as for the single-channel case, 
Eq. (2.5), but has to be constructed numerically [51-54]. It now depends on the constant 
normalization of the corresponding pion and kaon form factors T™ , s K (s = 0). In contrast 
to the vector case, the normalizations of the scalar form factors are not fixed by symme- 
tries; they are however related to the corresponding masses by the Feynman-Hellmann 
theorem [55, 56], 

and similar for the scalar form factors of the kaon. At leading order in the quark mass 
expansion, one therefore has r™(0) = 1, T*(0) = 0, F^(0) = 1/2, and 1^(0) = 1. Beyond 
that, information on these quantities can again be deduced from lattice calculations. We 
vary the normalizations of the kaon form factors according to T^(0) = 0.4 . . . 0.6, T S K (0) = 
0.95. . . 1.15, as suggested by the uncertainties in the corresponding low-energy constants 
given in Ref. [46], while keeping the rather well-known pion- form- factor normalizations 
fixed at T™(0) = 0.98, (0) = 0. 3 

In Figs. 3(b) and (c) we show the resulting invariant mass distributions for the pion 
pair for (A| = lGeV" 2 , A| = 0) and (X q s = 0, X s s = lGeV" 2 ), respectively, with the 
uncertainty bands as dictated by the above estimates for the uncertainty in the kaon form 
factor normalizations. In panel (6) the /o(500) (or a meson) shows up as a broad bump 
with a clear non-Breit-Wigner shape, while the /o(980) produces a peak exclusively in 
the strangeness form factor, panel (c). Thus, were a pronounced peak just below 1 GeV 
observed in r — > fj>ir + ir~ , it would allow one to straightforwardly extract X s s from data, 
without the need to employ any assumption on the internal structure of the /o(980) — 
additional information can be gained from also studying the KK final state, which however 
we will not detail here. This highlights the advantage of our approach compared to the 
one of Ref. [21]: in that work assumptions on the quark content of the /o(980) need 
to be employed in order to derive bounds, which then in turn strongly depend on these 
assumptions. In our case the bounds can be deduced directly from a fit to the spectra, 
once they are measured. 

3 A recent analysis of these scalar form factor normalizations directly based on lattice data, following 
the generalized framework of resummed chiral perturbation theory as in Ref. [57], yields F™(0) = 1.000 ± 
0.005, (0) = 0.013 ± 0.009, r£.(0) = 0.56 ± 0.06, and T S K {0) = 1.19 ± 0.11, thus perfectly compatible 
with the values assumed above, even though still more precise in some cases. We are very grateful to 
Veronique Bernard and Sebastien Descotes-Genon for communicating these results to us prior to publication. 
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4 Discussion 



If super symmetry was to show up in experiments like r — > /X7r + 7r~, there is no a priori 
reason why pion pairs in the vector channel would be significantly more populated than 
pion pairs in the scalar channel. While the effective couplings for the vector channel are 
given by squark exchange, the scalar channel is driven by sneutrino-exchange contributions. 
Thus one should expect interferences of the three currents discussed individually above. In 
this context it is important to stress that, below the first significant inelastic threshold, the 
phase of the form factor agrees with that of the elastic scattering amplitude [36], which is 
well known in both the scalar and the vector channel [39-41]. 

The Belle collaboration has given upper limits on branching ratios B(t~ — > n~Tr + ir~) 
with different kinematical cuts [58-60]. In particular, they find 



B(t~ fi-p° (770)) x B(p°(770) vr+vr") < 1.2 x 10" 
B(t~ /i~/ (980)) x 8(/ (980) -> tt+tT) < 3.4 x 10" 

B(t~ -> ^vt+tt") < 2.1 x 10" 



(4.1) 



The resonance signals are isolated by applying cuts to the tt + tt invariant mass spectrum, 
specifically 906 MeV < ^ < 1065 MeV for the / (980) [58] and 587 MeV < y/s < 962 MeV 
for the p(770) [59]. As we aim for deriving upper limits on coupling constants from null 
experiments, what will effectively enter the bounds on the scalar couplings are the lower 
limits of the uncertainty bands on the scalar form factors. Note that the vector form factor 
is known to a precision that any uncertainty therein is totally irrelevant at the accuracy 
we aim for. When comparing to the last, inclusive, bound in Eq. (4.1), we will set all 
form factors to zero above 1.1 GeV where we deem our representations not very reliable 
anymore; obviously, the bounds could be improved upon if lower bounds were available 
also at higher energies. Integrating Eqs. (2.9) and (3.4) in the respective ranges, we find 



/ 



962 MeV 



587 MeV 



( dyfs ' dy/s 



oTv dT s 



d\fs 



1065 MeV 



dTx 



+ {1.2... 1.8}|A^| 2 

< 2.7 x 10~ 20 GeV 
dT s 



0.94|Ay| 2 + {2.3. 
x 10" 3 GeV 5 



3.0}|Ag| 2 



{2.5...3.7}Re(AgA|*) 



+ 



<Vs= 0.10|Ay| 2 + {0.3. 



906 MeV \dy/s dyfs 

+ {0.35 . . . 0.42} Im(AgA!*) + {2.5 . . . 3.7}|A| 
< 7.7 x 10~ 20 GeV , 



0.6}|Ag| 2 + {0.5...2.1}Re(AgA|*) 



s |2 



x lO^GeV 1 - 



1100 MeV 



dV v dT s 



2Mt, 



dyfs= 1.05|Ay| 2 + {4.3. 



5.1}|Ag| 2 



ydy/s dy/s 

+ {0.47 . . . 0.57} Im(AgA|*) + {3.2 . . . 4.7}|A 



{2.4... 4.7} Re (Ag A|*) 



s |2 



x lO-^GeV 5 



< 4.8 x 10- 20 GeV 



(4.2) 
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product of couplings bound susy mass eff. coupling 



A 21i A 31i 


2.1 • 


10" 4 


m,- 

di 


Ay 


A 2il A 3il 


2.1 ■ 


io- 4 




Ay 


hi2^ai> A 2i3A-i! 


1.3 ■ 


10~ 4 


m Vi 


A§ 


A3i2A-2 2 , A 2 i3A-22 


1.5 ■ 


10" 4 


me. 


A s 5 



Table 1. Our best bounds on the products of supcrsymmetric R-parity- violating couplings. We 
assume all other couplings to be zero. The bounds scale with the square of the given mass divided 
by 100 GeV, so for a 1 TeV mass they are a factor 100 weaker. To obtain the bound wc employed 
the effective coupling in the far right column. Note that in the first two rows for the case of i = 1 
both operator products contribute in Ay and the relative phases are essential. They can only be 
treated separately if mj <C or <C . 

where the limits on the branching fractions have been converted into partial widths with 
the help of the lifetime of the r lepton. The uncertainties shown for the coefficients of the 
scalar couplings are due to the ranges assumed for the kaon scalar form factor normaliza- 
tions as discussed above. They are displayed explicitly in order to indicate the remaining 
potential for improvement, once this specific hadronic input is still better known. Note 
that contributions oc Im(AgA|.*) have to come from above the KK threshold, where the 
two interfering scalar form factors are not required to have identical phases according to 
Watson's theorem [36] any more. 

In order to set limits on the underlying coupling constants, we shall make the usual 
simplifying assumption that they are all real. Thus, assuming the first operator in Eq. (2.7) 
dominates we obtain for example the bound 

Too^v) ■ (4 ' 3) 

The corresponding bound on the second operator is 

A complete list is given in Table 1. Note that while the most restrictive bound on the vec- 
tor coupling Ay and the underlying products of fundamental coupling constants is derived 
from the dedicated search for r~ — > n~p° (770) [59], the best limits on the two effective 
scalar couplings stem from the most recent limit on r _ — > /i~7r + 7r~ without further kine- 
matical cuts [60] . In the literature the best previous bound is given from different processes 
bounding the individual couplings separately [16]. Combining them to a product bound 
we obtain 

where we have set the mass of the different virtual supersymmetric scalars equal. We have 
thus improved this bound by more than a factor of 30. We can also compare our result to 
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the related bounds obtained in Ref. [61]. The authors consider the effective operators 



O j = ^-(flT j T)(g a T jq ^, (4.6) 

independent of their origin, for different Dirac structures Tj. These contribute for example 
to r — > /i7r + 7r~ . Setting C 3 = Air the authors obtain the bounds 

A > 2.6 TeV (scalar) , A > 12 TeV (vector) , (4.7) 

where 'scalar' and 'vector' denote the cases Tj = 1 and Tj = j a , respectively. Comparing 
the vector case to our bound we thus have 



*21i*31i = 

4m 2 A 2 ' 

hence 

> ~ \ 2 

i' v — i I \ ^ o k 1 n 3 



(4. 



which is more than an order of magnitude weaker than our bound, partly owing to weaker 
experimental bounds at the time. We point out that in Ref. [61] the scalar form factors 
were assumed to be constant, which is, as we have discussed, not a good approximation. 



5 Conclusions 

Supersymmetry has to-date not been observed. Thus well-motivated versions, such as R- 
parity violation, which however have been considered less conventional, are now also being 
investigated in more detail. Most of the strictest bounds on supersymmetry with R-parity 
violation arise from precision, low-energy observables, e.g. meson decays, or decays involv- 
ing mesons. These computations typically involve simple approximations of the relevant 
hadron physics, such as the current algebra approximation. In this paper we demonstrate 
that the bounds on the R-parity-violating couplings can be considerably improved when 
including dynamical aspects of hadron physics. We do this in a model-independent man- 
ner, such that it can easily be applied to test other fundamental theories of physics beyond 
the Standard Model. 

To be specific, we have focused on the decay — > /i ± 7r + 7r~. Here, the hadron 
physics aspects can be treated particularly rigorously, as the strong final-state interactions 
of the pion pair are described in terms of the vector and scalar form factors, which are 
either directly measured experimentally, or can be reconstructed using the rigorous, model- 
independent methods of dispersion theory. We then employ upper bounds found by the 
Belle collaboration on the branching ratios of the lepton-flavor- violating r decays, and ob- 
tain bounds on products of the R-parity-violating couplings. Due to the extra information 
we have included and the slightly improved experimental data we find bounds which are 
more than an order of magnitude stronger than previous ones. 
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A Derivation of effective operators 

In this Appendix we present as an example the computation of the effective Lagrangian 
given in Eq. (1.9). We show the results in four-component spinor notation, as this still 
is the most widely used convention across communities. However, we actually derived 
the results in two-component spinor notation [62] and then translated them into four- 
component notation using Appendix G of Ref. [62]. 

The superpotential in Eq. (1.2) contains the terms e a b\ ; i j k L^Q b jDf : . In component fields 
the corresponding Yukawa Lagrangian terms are given by, see for example Ref. [63], 

-2lqd = -Kjk(d* kR vfPLdj - d* kR e?P L Uj + d jL d k P L ^ 

- Uj L d k P L (,i + Vid k P L dj - ludkPLUj + h.c.J , (A.l) 

where we have assumed that the Yukawa coupling is real and we have denoted the charged 
lepton by £j. The charge conjugate field is defined as l\ = Clf . If we now combine the 
second term for i = 3 as well as the hermitian conjugate of the second term for i = 2, 
identify the third index, and assume the intermediate d k squark to be very heavy, we 
obtain the effective Lagrangian 

^cff = -^4^ ^ P LUi) (u k P R p c ) . (A.2) 

dj 

A Fierz reordering as well as employing the definition of the complex conjugate field then 
gives the form for the effective interaction in Eq. (1.9) 

^cff = \'^^{u k 1 a P L u i ){jl la P L T) . (A.3) 
2 m- 

B Comparison to the narrow-resonance approximation 

In this Appendix, we briefly compare our treatment of the r — > /i7r + 7r~ decay in terms of 
form factors to the narrow-resonance approximation. If, e.g., the p meson were a stable 
particle, instead of Eq. (2.2) we had to use, cf. Ref. [17, 18], 

M = {p(k)\r\Q)M r Va , ( P (k)\j a \0) = (p(k)\Uui a u-drd)\0) = ±=f p M p e(\Y * , 

(B.l) 
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where f p = 220 MeV denotes the p decay constant and e(A) is the p polarization vector. 
The expressions for the physical and the stable p are most easily compared on the level of 
the decay rates. Employing the identity 

£e(A); e (AV = -^ + ^ ; (B.2) 

A P 

one finds that switching from a stable p to an unstable one is obtained by the replacement 
i/ p 2 M p 2 ^-i(2vr) 3 y ds j d<S> 2 (q^ + ^-)( P „ + - p^f\F v {s)\ 2 

= ^2 / ds@ ( s ~ 4M 2 )^±£ \F v (s)\ 2 , (B.3) 



where q 2 = s, we use the definition of the phase space as given in Ref. [19], and (p n + — 
p n -) 2 = — 4p* + 2 in the center-of-mass frame. In order to understand under which cir- 
cumstances the full expression is well represented by the left hand side of Eq. (B.3), we 
may parametrize the pion vector form factor by its spectral function, assuming a constant 
width. Then we may write 

As a consequence of unitarity the spectral function is normalized, 

dsp(s) = l. (B.5) 



oc 



LA/2 



In experiments the p meson is typically identified by imposing cuts on the two-pion invariant 
mass. If for the sake of the argument here we assume that these cuts were sufficiently wide 
that the normalization integral is exhausted, we find, using the explicit expression for the 
p width 



Lp ~ 6vT M 2 



as well as the KSFR relation g 2 pw7T = M 2 /(2f 2 ) [64, 65], 

f P = V, , (B.7) 

which is a known connection e.g. in the hidden-local-symmetry approach [66]. With f n = 
92.2 MeV, this identity is fulfilled to about 20% accuracy. However, in any realistic situation 
the p meson in the final state can only be isolated by cuts on the tttt invariant mass, such 
that the p spectral function is not fully saturated, leading to uncontrolled inaccuracies in 
the extraction of the effective parameters. 
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